Large voltage from spin pumping in magnetic tunnel junctions 
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Abstract 

We studied the response of a ferromagnet-insulator-normal metal tunnel structure under an 
external oscillating radio frequency (R.F.) magnetic field. The D. C. voltage across the junction 
is calculated and is found not to decrease despite the high resistance of the junction; instead, it is 
of the order of //V to 100//V, much larger than the experimentally observed value (100 nano-V) 
in the "strong coupled" ohmic ferromagnet-normal metal bilayers. This is consistent with recent 
experimental results in tunnel structures, where the voltage is larger than //Vs. The damping and 
loss of an external RF field in this structure is calculated. 

PACS numbers: PACS numbers:73.40.-c,71.70.Ej,75.25.+z 
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There has been much recent interest in the spin dynamics in hybrid nanostructures com- 
posed of ferromagnetic and normal metal layers. [1-5] Following earlier work on the spin 
torque effect, the spin pumping effectfl has been de—tedfor "strongly coup.ed" oh m ic 
metallic multilayers as an additional contribution to the FMR linewidth in FM/NM multi- 
layers (where NM is Pt, Pd, Cu, etc.) [3] and more recently as a difference in voltages of the 
order of 100 nV between two FM/NM interfaces of a NM1/FM/NM2 trilayer.^] Two types 
of metallic structures are commonly studied. In addition to the strongly coupled ohmic mul- 
tilayer systems, weakly coupled tunnel structures have also been extensively studied. The 
physics of these two types of systems can be very different. 9[] In particular, for the tunnel 
structures, the coupling of the longitudinal magnetization and the charge lead to magne- 
tization and charge dipole layers at the interface. After the effect of the electron-electron 
interaction is included, it is found that because of the large difference of the length scales 
associated with the charge (screening length, lA) and the spin (spin diffusion length, 
~ 100 A) fluctuations, there is a larger splitting of the chemical potentials than that pre- 
dicted by the conventional spin accumulation picture. flo| Whereas the conventional picture 
suggests that the splitting scales with the current and will decrease with an increase with 
resistance, this is no longer true in the more complete picture. 

Recently Moriyama and coworkers ll|] reported measurements of the dc voltage attributed 
to the spin pumping effect in different tunnel junctions, and demonstrate that the voltage 
is larger than micro-volts, enhanced orders of magnitude compared to that for metallic 
trilayers. In this paper, we generalize our recent work on spin torque 13|] to the spin pumping 
situation and found an enhanced voltage for the tunnel structures, in agreement with the 
experimental results. We now describe our results in detail. 

The system we have in mind is a ferromagnet-normal metal tunnel junction where the two 
interfaces between the ferromagnet-insulator-metal sandwich structure are assumed to be at 
z = ±d/2. We assume the z axis to be perpendicular to the faces of the tunnel junction. 
The initial magnetization is assumed to be in the x-y plane with an orientation given by 
Po = e x f° r the ferromagnet on the left hand side of the sandwich structure. 

Because the work functions of the metals on opposite sides of the junction may not be 
equal, at zero external radio frequency (RF) field there will be a charge dipole layer formed 
at the interfaces. What we are calculating here are the changes from the zero field situation. 
This surface inhomogeneity can lead to an additional contribution to the increase in the 
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FMR damping, as we explain below. The experimental structures usually possess edge 
domains where the switching starts. The magnetization is thus not completely uniform in 
the x-y plane. To bring out the essential physics, we shall not consider this complication in 
the present paper but we hope to come back to this in the future. 

Under an external time varying RF field, we expect the magnetization in the ferromag- 
net to be a sum of a uniform magnetization which is a solution of the inhomogeneous Bloch 
(Landau-Gilbert) equation due to the external field and a spatially varying solution of the 
homogeneous equation so that the boundary condition can be satisfied. This spatially vary- 
ing part provides for the additional damping and the voltage observed in the experiments. 
Our approach is to obtain general solutions in each part of the junction (Eq. (11), (12), 
(15), (16)). The amplitudes of these solutions are determined by the boundary conditions 
(Eq. (7)). From these amplitudes, the voltage and the damping can be derived. We first 
describe the general solution of the magnetization in a ferromagnet. 

Magnetization in a ferromagnet: Our starting point is the equation of motion of the 
charge and the magnetization. For the charge, it is just the equation of charge current 
conservation 

r-r . 95n 

V-Je = -g r (1) 

where J e is the total current. The equation for the magnetization M has been much dis- 



cussion extensively in the past. 12| The equation takes the form of the phenomenlogical 

classical Landau-Lifshitz (Bloch) equation with longitudinal and transverse damping and an 

additional source term 

dM - <M 

— 7 M x H - aM x (M x H) + V • J M = (2) 

CjTj T 

where 7 is the gyromagnetic ratio, and H is the effective field describing the precession of 
the magnetic moments given by H = H e + H an + H dip + H cx . H cx = JV 2 M is the effective 
field due to direct exchange; the anisotropy term includes a bulk and a surface anisotropy 
energy H an = H ab + H as , H ab = KM , H as = K S M S where M s = M 5(z + d/2). Here we have 
separated a bulk and a surface contribution that acts on the surface magnetization M s . For 
simplicity we have assumed this surface contribution localized at the interface. The other 
terms can also contain a surface contribution and can be treated in a similar manner as this 
ansiotropy contribution. For simplicity of presentation we illustrate our results with just 
this term. H e represents the external field; and H dip denotes the dipole-dipole interaction. 



Jm is a spin current (tensor). The currents are driven by density gradients (diffusion) and 
external forces. 

J e = -aVV - eDVSn - D M V(£M • p ) 

Jm = -<t m V(Ypo) ~ £>mVAM - D'V(5np ) 

where a, om are the effective conductivities for the charge and magnetization. p is a unit 

vector along the direction of the equilibrium magnetization: p = °, with M the local 

|M | 

equilibrium magnetization, AM = M(l — |M|/M ) is a change in magnetization. D, D', 
D M , D' m are the effective diffusion constants. V = V e + W, with V e the electric potential 
describing the external electric field and W the local electric (screening) potential due to 
the other electric charges determined self-consistently by 

W{r) = J dh'Uiv - v')5n{r') (4) 

with U the Coulomb potential. The total number density of charge carriers and x component 
of magnetization are given by n = J2 S n s, M x = J2 S sn s- In the coordinate system with one 
of the coordinate axis along the direction of the magnetization, the spin current can be 
understood as the difference of the spin up current and the spin down current. The vector 
dependence is such that the equation is covariant. The Landau-Liftshitz equation without 
the source term V • Jm is believed to describe the physics of ordinary domain walls where the 
direction of the magnetization changes but its magnitude remains fixed. Eq. (3) is consistent 
with this belief. For ordinary domain walls, Jm = 0. r is the longitudinal relaxation time, 
describing the relaxation of the system towards its local equilibrium value of magnetization. 
a measures the transverse (Gilbert) damping term. 

Substituting the expression for Jm into the modified Landau-Lifshitz equation (2) we 
obtain the linearized relaxation equation for M: 

V 2 <M- (-^-+iu/D' M )M+{Po x (V 2 cM- ^-5M-K s 5M s /l 2 sf ) + 7 (<5M x H +M x H x ) (5) 

L sf L sf 

+aSM ± M H /D' M = -(£>' /D' M )p (V 2 6n - -^) 

where only H cx = JV 2 M = JV 2 cM and H an = K~M are kept in the precession term 

7M x H, and use has been made of Gauss' law: V 2 V = V 2 W = 8n. The bare spin 

diffusion length l s f and the bare screening length Ao are given by / 2 j = tD' m and Aq = 

respectively. Other dimensionless parameter are C = ^ and K = M. 

Dm J 



The charge current conservation (1) yields, for the steady state without linearization, 

-^5n - (V 2 - iu/D)5n - (D M /D)V 2 {6M ■ p ) = 0, (6) 

which, together with Eq.(jS]), describes the distribution of the charge and magnetization away 
from the tunnel junction in terms of their values at the junction. (To simplify the algebra, 
we have made the approximation that D=D'). The values of the charge and magnetization 
densities at the junction can be determined by matching boundary conditions across the 
barrier. We first solve these equations in the metal part of the junction. These solutions 
determine the charge and magnetization dipole layers. 

The solution of eq. ([5]) can be written a sum of two terms, 

5M = 5M + 6m\ 

a spatially uniform (5m l ) solution of the bulk inhomogeneous equation with the source 
term M x and a sptially varying solution (5M ) of the homogeneous equation. The 
inhomogeneous bulk equation is 

1 K 
-{ w + (iu - a')/D' M )5m - Cpo x + ^( 6m x H + M x Hi) = 

L sf L sf 

where a' = aM^Ho. This is the conventional FMR equation, which can be readily solved. 
Define e± = e z ± ie y , then x e± = ±ie±. We write the transverse magnetization as 
5mY = J2 Sm±e± and obtain 

Sm± = x±Hi,± 

1 K 

where l/x± — [( — ±i-f2~ + (iu} — a')/D' M ) + (-^- + jH ]/'yMQ. Associated with this transverse 

\sf l sf 

magnetization, there is a change of the longitudinal magnetization given by 

5ml = M - (Mq - 5m l l) 1/2 w 0.55m^/M . 

This is the lowest order correction to the longitudinal magnetization. Higher order non- 
linear corrections to the transverse magnetzation will produce changes in the longitudinal 
component that is higher than 3rd order in H\. In the equation of motion (2), no lower order 
correction are produced. 

The equation for the spatially varying term becomes 

V 2 M - (i + iu/D' M )6Nl + Cpo x [V 2 cM - -^<M - k s (6ui* + 5M 0s )/l 2 sf ] + 7 5M x H 
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+a5M 0± M H /D' M = ~(D' /D' M )p (V 2 5n - 



The solution of this equation is similar to that in our previous studies. 13|] Away from the 
boundary, the surface terms are absent. This equaion becomes homogeneous. We solve 
this homogeneous equation and form linear combinations of them to satisfy the boundary 
constraints. By integrating this equation over a small region of space at the boundary 
we arrive at the condition that the difference between the tunnelling and the ferromagnet 
pseudo spin current is equal to surface anisotropy term: 

JL-Im = JK S M x (cW + 5M 0s ). (7) 

where is the tunnelling magnetization current, the pseudo spin current [jjj] Jm = Jm — 
7 JMo x d z 51SA includes an extra term involving the exchange that affects only the transverse 
magnetization current. We expect this extra term to be also present for ohmic junctions 
but so far it has not been included. In previous spin pumping studies on ohmic junctions, 
a term of a similar functional form gn x dn/dt (n = M/|M|) has been discussed. However, 
the coefficient was interpreted as a spin mixing conductance. We next discuss the solution 
of the homogeneous equation. 

We expect the charge and magnetization dipole layers to decay away from the interface 
with length scales controlled by the spin diffusion length and the screening length. Because 
of the vector nature of the magnetization, there are three normal modes by which they can 
decay away from the interface. Including the charge degree of freedom, there are four normal 
modes that one can consider. For the ferromagnetic metal on the left hand side, we thus 
consider the following ansatz: 

4 4 



6n L = J2Snf e {z+2)/l \ <Mq = £<M&e ( * +S)/,< , (8) 



i=i i=i 
where the superscript L denotes the left hand side. 

Letting the coefficients before the exponential scaling functions vanish for steady-state 
solutions, we get for small uj the renormalized screening length 

h = Ao£i /2 , (9) 

the renormalized spin diffusion length 



and a combination of the exchange length and the spin diffusion length 

h,A = " ±*C«r)/(l " ±<) + (iw - a')l 2 sf /(D' 2 4 (l - ±<))] 1/2 . 

The £s and (3 are measures of the asymmetry of the spin up and spin down conductivities 
of the ferromagnet: & = [1 - D' D M /{DD' M )}/[1 - a M D'D M /(aDD' M ) + iuj\ 2 /D}. £ 2 = 
(1 - /3 2 )/(l - iul 2 sf /D' M ), f3 2 =[l- D'D M a M /(aDD' M )}. k t = k + 1 H l 2 sf /(M . As we shall 
see below, I3 and U correspond to length scales with which the "precession" dies away from 
the interface. The additional term / y5M x Hq modifies these two lengths accordingly. The 
screening length and the spin diffusion length are renormalized. From eq. (6) we find that 
the charge densities can be related to the magnetization densities by 

6n L l0 = e(ef - K = J^E m 20, *»» = *4 = 0. (10) 

Because h » Ao, 5ri2o/e « 5M 2 o/fiB- As we see below, in general 5M 20 is much less 
than 8M W . Inserting the "eigen-solutions" into equations (jHJ), we finally obtain analytic 
expressions for the dipole layers: 

Sn L = 6n^ z+ ^ h + 5n^ z+ i )lh (11) 
SM L = p%5Mf e iz+ ^ /h + p^5M^ z+ ^ h + e^M^e^^ + eHM^e {z+ i )/h . (12) 

The two transverse modes corresponds to the left and right circularly polarized modes e±. 
8M?q, with i = 1,2,3,4, are to be determined later. Terms of the order (X /l s f) 2 or higher 
have been neglected since l 2 j » Aq- Also, to simplify the algebra we have assumed that the 
ferromagnetic thickness dp to be larger than the spin diffusion length so that we do not need 
to worry about "reflection" effects from the leads. As advertised, the charge dipole layer is 
the sum of two terms, one decaying with a length scale of the screening length; the other, 
the spin diffusion length. The vector magnetization dipole is now a sum of four terms. The 
first two ( (SVI^q, (SVIgQ ) are along the direction of the original magnetization; the last two 
are perpendicular to the direction of the original magnetization and describes the precession 
of the magnetization around the original axis. Again, the first two terms correspond to 
decay lengths of the order of the spin diffusion length and the screening length, while the 
precession term only decays with a length scale that is a combination of the exchange length 
and the spin diffusion length. 
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With equations (iTTj) and (|T2j) . the charge and magnetization currents J e and Jm can be 
worked out as 



Je = ^E ext (13) 

Jm = CTjwE ea; tPo /e + 

+ (hM^L^ + ^k e ^ e+5M L e ^3 + Pk eze _ 6M L Qe 5 Z / U (14) 

'2 '3 '4 

where 5z = z + ~, E e;rt = E ext e z is the external electric field inside the conductor. Note 
that the magnetization current is not a function of the rapidly varying part of the charge 
and magnetization densities Snio, 8M\q. In principle, the magnetization current can contain 
a term of the form J mi exp(z/li). In the generalized Landau-Gilbert equation (Eq. (2)), 
terms of different functional dependence are each equal to zero. The only terms that are 
proportional to exp(z/Zi) comes from V • Jm and is proportional to Jmi/1\- This term and 
hence its contribution to the magnetization current is equal to zero. To match the quantities 
at the boundaries we next consider the charge and magnetization in a normal metal (N). 

Normal metal: On the N side, the charge and magnetization are not coupled. The charge 
is given by 

5n R = <5n^exp(— z/X). 

The magnetization satisfies the equation c^M = {D n d 2 — l/rj^)M = 0. From this we obtain 

SM R = 5M R exp[-(z - d/2)/l R \. 

The longitudinal magnetization current at the interface (z=d/2) is given by J M = 
-D N 5M*/l R . 

The longitudinal magnetization current at the left interface is given by eq. (16). Equating 



J M to J M , we get 



5M* = -(1 - (3 2 )D M 5M£ l R /[D N l 2 }. (15) 



The magnetization on the right is proportional to SM^q and is not a function of SM[ . As 
we shall see below, SM^ Q » SM^,, hence the longitudinal magnetization change on the 
right is much less than that on the left at the boundary. The charge neutrality condition 

_d 

jd° 6n L dz + 5n L dz = yields 

2 

8n*=-(l 1 8n 1 S ; + l 2 5n% )/X, (16) 



S 



These two equations express the quantities on the right in terms of quantities on the left. 
We now determine the amplitudes of these physical quantities by matching the boundary 
condition as in eq. (7). 

Boundary conditions: The longitudinal magnetization current in the ferromagnet arriving 
at the interface I^is equal to the magnetization current J l M across the interfacce due to 
tunnelling because the term on the right hand side of eq. (7) is along the transverse direction. 
For the longitudinal component, I m — J m - The longitudinal magnitization tunnelling current 
is equal to the difference of the spin up and the spin down tunnelling current. From standard 

n 

calculations of the tunnelling current [151] we get 

Jm = E s\T ss ,\ 2 {8n Ls {E + Srf) - 5n Rs (E + Sfif)). 

s 

Here S/j, contains contributions from the electric potential due to the charges at the interface 
and that from the accumulation due to the bottleneck effect. The change of the electron den- 
sity of spin s can be related to the change of the total charge and magnetization densities by 
(we use units so that \is = 1): 8n s = 0.5(Sn+s6M x ). The longitudinal magnetization density 
is the sum of contributions from the solutions of the homogeneous and the inhomogeneous 
equations: 

5M X = 5mi + 5M x0 . 

From eq. (14) J\ L = (1 ~^ )O W M 2 l . The inhom geneous term 5m t is uniform and does not 
contribute to the magnetization current inside the ferromagnet. From J l M = Jf f , we 
get 

(1 - (3 2 )D^5M 20 /l 2 = s\T ss \ 2 (6n Ls - n Rs ). 

s 

All variables of this equation can be written in terms of the two independent variables 
<5Mi 0i 2 - Now 5n 0Ls = 0.5[5n l0iL + 5n 20 ,L + s(5M l0iL + 5M 20 ,l)-] Using eq. (11) and (16) we 
get 

(1 - (3 2 )D'^5M 20 /l 2 = ]T s\T s8 \ 2 [8M^{e - 1 + s + h(£ L - 1)/A) + sSrri]. (17) 

s 

This equation implies that 8M 20 is of the order of c t SM^ /c m where c t (c m ) is the tunnelling 
(metallic) conductance. c t much smaller than the metal conductance c m . Thus SM 2Q is much 
smaller than 8M±o. 

For an open circuit, the total charge tunnelling current is zero. We get J = 
J2 S \T ss \ 2 (8nL S — riRs) = 0. Substituting in the expresssions for the charge densities and 
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using the condition that <5M 20 << 5M W , we get 

£ \T ss \ 2 [8M^{i L -1 + 8 + h(£ L - 1)/A) + s<5t<] = 0. 

s 

Solving this equation, we finally obtain 

5M 1 L = -fSmi (18) 

where / = (£ s \T ss \ 2 s)/[E s \T SS \ 2 {£ L - 1 + s + h{£ L - 1)/A)]. The corresponding charge 
is, from eq. (10), 6n^ = (£ L — l)/5m*. The charge and the magnetization densities are 
proportional only to the ratio of the conductances. Hence they are not necessarily small 
for tunnel junctions. As we emphasized before [9j], this comes about because A << l a f. 
Emf: 

The DC voltage is estimated as the change of the mean chemical potential across the 
interface, given by AV = 0.5 X) s A(<5/i s ) = 0.5 £) s A(5n s /iV s ) where N s is the density of 
states. This drop includes a contribution from an electric potential as well as a contributions 
from electron density changes due to bottleneck and electron-electron interaction effects. 
This drop can be written as (5M L » 5M R ) 



AV = 0.25(e5mi/fi B )f[^ L - 1)(1/JV* + 21 1 /{X R N R ) + 1/iVf) + 1/N^ - 1/JVf]. 

The longitudinal magnetization density is 5m l x = 0.5(5m t 1 _) 2 /(Mqv ) = 0.56 2 Mq/v where 
v is the atomic volume, 9 = 5m 1 / M is the precession angle. Hence 

AV = 0.l25e6 2 M /(vfi B )f[(t L - 1)(1/N^ + 2h/(XN R ) + l/N L ) + 1/N$ - 1/JVf]. (19) 

As expected, this d.c. voltage is proportional to 8 2 , as is observed experimentally. Most 
importantly, it is proportional only to a ratio of the conductances. Hence its magnitude 
is not small. The factor f, as given after eq. (18), depends on the asymmetry between the 
majority spin and the minority spin conductances in the insulator. The larger the difference, 
the larger the value of |/|. We next estimate the order of magnitude of DeltaV. 

We expect Mq/(vhb) to be of the order unity, e/N ( N is the average density of states ) 
to be of the order of 0.1 volt. Dependening on the asymmetry between the majority and the 
minority spin tunnel conductances in the insulator, the value of f can range between 1 and 
0.1. Similarly, depending on the asymmetry between the majority and the minority spin 
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conductances in the ferromagnet £ L — 1 can range in value between 1 and 0.1; 1/N+ — 1/AL 
to be of the order of 1/N to 0.1/N. Hence AV « (10~ 2 - 1O _4 )0 2 uoZ*. For « 0.1, « 
(10~ 4 — 10~ 6 )t> o/t, in agreement with the experimental results, which is larger than microvolts. 
We next address the issue of damping. 

Damping: The loss can come from three sources: (1) from the interface inhomogeneity, 
(2) from loss of the transverse magnetization current through the barrier, (3) from loss of the 
longitudinal magnetization current. As we explain below, these contributions have different 
dependence on the external RF magnetic field. The contributions for the first two sources 
to the damping coefficient are independent of the field strength; that from the last source is 
proportional to the input power. The contributions from the last two sources are inversely 
proportional to the junction resistance and thus are much smaller for tunnel junctions. 

We first estimate the loss connected with the longitudinal magnetization. This loss is 
equal to J2 S i1 r 's where j s , r s is the current and junction resistance for spin s. This is of the 
order of (Sm l x ) 2 \T\ 2 . Since 5m x is proportional to the input power, this loss is proportional 
to the power squared. Its contribution to the damping coefficient is obtained by normalizing 
the loss by the energy density and hence is proportional to the power. Because this loss 
is proportional to |T| 2 its contribution is small for tunnel junctions. Similarly, we expect 
the transverse magnetization current to incur a loss of the order of {5m\) 2 \T\ 2 . Since Sm±_ 
is proportional to the field, this loss is propotional to the power. Its contribution to the 
damping coefficient, again obtained by normalizing with espect to the energy density, is thus 
independent of the power. This loss is also proportional to \T\ 2 and will be small for tunnel 
junctions. 

We next estimate the loss connected with the interface inhomogenity. This requires 
knowledge of 5M3o,4o which we now determine. Again, we expect the transverse magnetiza- 
tion to be a sum of a term that is the solution of the inhomogeneous equation (5m 1 ) and 
terms that are solutions of the homogeneous equation (<5M 3j4 ). We calculate <5M 3j4 using the 
boundary condition given by eq. (7). From eq. (15) the transverse magnetization current 
at the boundary is 

d' n' 

J M — 7 e + OM 30 ; e_OM 40 

i 3 t 4 

The pseudo spin current in eq. (7) is thus given by 

Im = -{D' M + ilJM )5M^e + /h - (D' M - i 7 JM )e_5M 4 L /Z 4 
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Eq. (7) also involves the tunnelling transverse current J l m . To evaluate this we follow 



standard practice 



151 ] and calculate the rate of change of the transverse magnetization due 



to tunnelling. We found that the tunnelling current for the transverse magnetization can be 
written as J\ I+ = M+(gi + ig 2 ) + M+(g 3 + ig^j where 1(| #1,2,3,4 are proportional to |T| 2 . 
A similar equation for Jm- can be written down. This shows that the contribution from 
the tunnelling current is smaller than the other terms in eq.(7) and thus will be treated by 
perturbation theory. We finally obtain to lowest order — = 7fC s M x (Sm l + 5M 0s ). 
Substituing in the expression for I m , this equation becomes 

[±zD' M /(l±'yK a M Q ) - 1 - J/(l ± K a )]8M^ = 8m^. 

Here l + = I3, /_ = U, 5M + = 5M 3 o, 5M_ = 5M^. As we go away from the interface, the 
transverse magnetization density dies off exponentially. The total magnetization is given by 
5M±l± = 5m±/X± where X± = [±iD' M / \ , yK s M ) —l± — J/K s ]. The correction term due to 
the tunnelling magnetization current is equal to — X^I^/ (jK s M ). 

The magnetic susceptibility, given by (SM±l/dp + 5m±) / Hi t ±, becomes \ — X°(l + 
X^ 1 /dp). The additional damping comes from the imaginary part of x which now contains 
a term proportional to Re(x°)lTnX~ 1 /dp- This term is proportional to the metallic "resis- 
tance" D' M which in turn comes from the spatially varying part of the magnetization induced 
by the surface, as we have anticipated. This contribution is not a function of the junction 
resistance and will be of the same order of magnitude for multilayers as well as for tunnel 
barriers. 

In conclusion, we discussed in this paper the voltage and the damping of an RF field 
in ferromagnetic tunnel junctions. The voltage is controlled by changes of the longitudinal 
magnetization whereas the damping seems mainly associated with the transverse magne- 
tization. Additional sources that can induce transverse magnetization localized near the 
interface can come from localized changes of the Hamiltonian such as the surface anisotropy. 
The calculation in this paper can be trivially extended to junctions with ferromagnets on 
both sides. For junctions involving two ferromagnets on opposite sides (F1-I-F2 or F1-F2), 
the interface anisotropy K s will contain a term from the dipolar interaction between Fl and 
F2. The loss will then be a function of the orientation of the magnetizations of Fl and F2, 
consistent with experimental results. 
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